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Abstract 

Using partonic states we verify the transverse momentum dependent factorization for single trans- 
verse spin asymmetries in semi inclusive DIS at leading but nontrivial order of a^. The factorization 
has been previously derived in a formal way by using diagram expansion at hadron level. We find 
with our partonic results that two relevant structure functions satisfy the factorization. Our results 
also satisfy the collinear factorization but with the perturbative coefficient different than that derived 
formally. 



1. Single transverse-spin asymmetries(SSA) have been observed in experiment. An updated review about 
the phenomenology of SSA can be found in [1]. These asymmetries are expected if scattering amplitudes 
have nonzero absorptive parts. SSA in a scattering of elementary particles is easy to be understood, 
but SSA in a process involved with hadrons, especially, with a transversely polarized hadron, is difficult 
to be predicted, because we have not enough information about the structure of hadrons. However, 
the observed asymmetries provide a new tool to study the structure because they are sensitive to the 
correlations of partons inside a hadron and the angular momenta of these partons. 

In general it is difficult to predict the size of SSA except some cases like production of a heavy quark 
with transverse polarization, where SSA can be studied by using perturbative QCD directly [2]. In the 
cases studied in experiment with an initial hadron being transversely polarized, it is not possible to use 
perturbative QCD directly. One needs to rely on QCD factorization in which perturbative- and nonper- 
turbative effects are separated. The nonperturbative effects are represented by various matrix elements 
of QCD operators and these matrix elements contain information about the structure of hadrons. Two 
factorization approaches have been proposed. One is by using transverse- momentum-dependent (TMD) 
factorization, where one takes transverse momenta of partons in hadrons into account. Another is the 
collinear factorization. It should be noted that in the two approaches the factorization is derived or 
proposed rather formally in the sense that one works at hadron level by using the diagram expansion. 
The two approaches are applicable in different kinematic regions. The TMD factorization can only be 
used for the kinematic region where the observed transverse momentum is much more smaller than some 
large scales, while the collinear factorization can be used if the observed transverse momentum is much 
more larger than Aqcd- The purpose of our work is to examine the two factorization approaches for SSA 
in Semi Inclusive DIS(SIDIS) at parton level, following our work for SSA in Drell-Yan processes[3]. 

In the approach of TMD factorization, the nonperturbative effects responsible for SSA are repre- 
sented by matrix elements containing T-odd- and spin- flip effects. These matrix elements are the so- 
called Sivers[4]- and Collins[5] functions. So far TMD factorization has been examined carefully only 
for physical quantities which do not contain T-odd effects [U [71 [HI O [10]. TMD parton distributions 
entering the factorization can be defined with QCD operators consistently. Intensive efforts has been 
spent to study how to consistently define or interpret Sivers function as a parton distribution which is 



gauge invariant and contains initial- or final state interactions [5} [TT\ \12 \ 1131 114j . Through these studies 
consistent definitions of these distributions or fragmentation functions which contain T-odd effects can 
be given. In [TOj the problem of universality of T-odd parton distributions and T-odd fragmentation 
functions has been studied. The TMD factorization of Sivers effect has been examined in a spectator 
modelfH]. SSA has been studied extensively in terms of Sivers functions [El [lU [171 [IBl [l9]. In the 
approach of collinear factorization, the nonperturbative effects responsible for SSA are represented by 
twist-3 matrix elements|20 t I2 H 122 ^ I23j. or called ETQS matrix elements. These twist-3 matrix elements 
contain only the effect of spin- flip interactions. The nonzero absorptive part or T-odd effect is generated 
by poles of parton propagators in hard scattering. Applications of the collinear factorization for SSA can 
be found in \24\ [25] . It is interesting to note that for SIDIS and Drell-Yan both approaches are applicable 
for certain kinematic region and it can be shown that the two approaches are equivalent [261 [27] . 

As mentioned, the above factorized results are derived rather formally in the sense that one works at 
hadron level by introducing various parton density matrices of hadrons in a diagram expansion. It should 
be noted that QCD factorizations, if they are proven, are general properties of QCD Green functions. It 
means that the two factorization approaches, if they hold, they should also hold by replacing hadrons 
with partons. It is the purpose of the study presented here to show how SSA in SIDIS can be factorized 
in the TMD-factorization approach by replacing hadrons with partons. In our previous work[3J we have 
shown that the TMD-factorization for SSA in Drell-Yan processes holds with partonic states at leading 
but nontrivial order of a^. With the partonic states one can also derive a collinear factorization for SSA 
in Drell-Yan processes but with perturbative coefficients different than those derived before. For SIDIS 
we replace the hadron in the initial- and final state with a quark q. In order to have spin-flip we keep the 
quark mass as nonzero and every quantity is calculated at leading power of m. It should be emphasized 
that the derived perturbative coefficients in the factorization formulas do not depend on the quark mass 
m. Two structure functions for SSA in SIDIS are studied in this work. 



2. We consider the SIDIS: 

e + hA{P, S^)^£ + 7*(g) + hA{P, S^)^i + X + hB{Ph)- (1) 

We take a coordinate system in which moves in the z-direction and the virtual photon moves in 
the — z-direction. The initial hadron hA_ is transversely polarized. We will use the light-cone coordinate 
system, in which a vector o'^ is expressed as a'^ = (o"*", a~ ,a±) = {{a^ + a^)/y/2, {aP — a^)/\/2, , a^) and 
a\ = {a^Y + («^)^- We also introduce two light-cone vectors: n'^ = (0,1,0,0) and l^^ = (1,0,0,0). The 
relevant kinematic variables can be defined as: 

= 77^ = 77^ . Zh = — . (2) 

IF ■ q IF ■ q F ■ q 

The hadronic tensor for the process is defined with the electromagnetic current J'^ of quarks as 

^'^^ = 4i;;E/ -0^J'-''{hA{p.s^)\j^{x)\xPh){xPh\rmhA{p,s^)). (s) 

We will be interested in the kinematical region of P^j^ ^ Q^, The hadronic tensor in this region has the 
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following leading- twist structures: 



(1) 



I ai^ pM 

(3) 
T 
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n(2) 



(4) 



with the notation g'^" = gi"" - nH"" - n"" l^" and e^^ = e"''^'"'lanp. Ph± is the unit vector Ph± = Ph±/\Ph±\- 



All structure functions F's depend on xb, z^, P?<, and Q^. The differential cross section is given by: 



da 



dx B dydzh d(j)s d'^ Pf 



hi. 



F^j^^ + sin( 



+(1 - y)xB - cos{2^h)Flj^ + \s±\ sm{^h + ^s)F^' 



p{2) 



+ -|sj_|sin(3(/)/i 



QF^ 



(5) 



where 5 is the invariant mass of the initial lepton and the initial hadron, y is the fraction of the lepton 
energy loss, (ps and (ph are azimuthal angles for the transverse polarization vector of the initial hadron 
and the transverse momentum of the final state hadron, respectively. A more general decomposition of 
the above differential cross section can be found in [28]. By studying the angular dependence of the 
differential cross section in experiment, one can measure the different structure functions. In this letter 
we will focus on F^^ and f!^\ 

It has been suggested that the two structure functions can be factorized with various TMD partons 
distributions and TMD fragmentation functions. The relevant TMD parton distributions and fragmen- 
tation functions can be defined by introducing a gauge link along the direction = {u^ ,0,0): 



L„(iboo, z) = Pexp ^—igs J dXu ■ G{\u z)^ . 



(6) 



Two relevant parton distributions can be defined in the limit -^u [51 [HI [13] : 

1 (' dz~d^zx_ 



q±{x,k±)e'^s±f,k±^ 
SqT{x,k±)s'^ 



2 J (27r)3 
1 f dz^d^zx_ 
2 



-ik-z 



-ik-z 



{P, s^mz)LU^, z)j+Lu{oo, 0)^PmP, 



(27r)2 



{P, s^mz)Li{oo, z)7+7^75L„(oo, 0)V(0)|P, s^), (7) 



with z^ = (0, z~ , Zj_). X is defined as = xP^ . In the above one only takes the spin-dependent 
part of the matrix elements into account. q±{x,k±) is the Sivers function[l]. A relevant TMD parton 
fragmentation function is defined with the gauge link L^(— oo,0) along the direction v'^ = {v^ ,v^ ,0,0) 
with <^ v^: 



{x,Ph±)pr 



hi. 



1 

4x 



dz^d^z_i_ 
~(2^ 



exp 



-iz 



-Ph/^) 



5^-Tr Uil^^T) {0\L,{-^,0)i^{0)\PhX){XPh\i>{z)Ll{-<x,,z 



X 



(8) 



3 



with = (z"*", 0, The parton carries the momentum k^^ = (0, /x, 0, 0) and the hadron is with the 
momentum Pj^ = (Pfj^ , Pf^ , Ph±) with as the large component. The function 5q is called as Collins 
function[5]. It should be noted that one can also define the Collins function in a frame through a Lorentz 
transformation where the hadron has zero transverse momentum and the parton has the transverse 
momentum —xPfi± [6]. The Sivers- and Collins function can only be nonzero if a nonzero absorptive part 
exists in the corresponding hadronic matrix elements. The definitions of other TMD parton distributions 
or fragmentation function of a unpolarized hadron can be found in [6l|8l|9j. One should keep in mind that 
all of these TMD parton distributions and fragmentation functions are defined in non-singular gauge. In 
principle in these definitions one should add a gauge link in the transverse direction at the space-time 
boundary as shown in [12], e.g., in Eq.(8) one should add a gauge link at x~ = oo along the transverse 
direction to make the definition gauge-invariant explicitly. However, in a non-singular gauge like Feynman 
gauge the gauge fields at infinite or x~ are zero. Hence the gauge links at infinite x"*" or x~ are unit 
matrices. In a singular gauge like light-cone gauges these gauge links can not be omitted [12]. We will 
work with Feynman gauge. The structure functions can be factorized as: 

A±) q±ixB,k±)q{zh,p_L)6'^{zhk_L +p± - Ph±), 

Phi.^ SqTixB,k±)Sq{zh,p_L)S^izhk± + p± - Ph±), (9) 

where q is the standard TMD parton fragmentation. Beyond the leading order of one has to implement 
a soft factor representing effects of soft-gluon radiation. 

For our purpose we replace each hadron state with a one-quark state. We consider the SIDIS of the 
quark states as: 

q{P,s^)+j*{q)^q{PH)+X, (10) 

with P^ and Pf^ as large momentum components. In order to generate SSA and nonzero q±, there 
must be exchange of two gluons at the leading order. One of the two gluon must be in the unobserved 
state X because P^ has nonzero transverse components. The exchange of another gluon must introduce 
a virtual contribution and generate a nonzero absorptive part in the scattering amplitude. With this 

(1 2) 

requirement one easily finds nonzero contributions to Fj. ' come from the interference of amplitudes 
given by diagrams in Fig.l. 

In Fig.l the diagrams in the first raw are the scattering amplitudes which can have a nonzero ab- 
sorptive part. The absorptive parts are determined by physical cuts which are given in these diagrams. 
To calculate the relevant hadronic tensor in the full kinematic region one has to calculate all diagrams 
in Fig.l. This will be very tedious. Fortunately we only need the leading contribution in the expansion 
of Ph±/Q'^- As shown in [3], one can perform the expansion before the loop-momentum integration of 
the two exchanged gluons. Before presenting our detailed results we point here that with our partonic 
state Fj. does not receive any nonzero contribution at the orders we consider. It is worth to point out 

('3') 

that the parton distribution function related to Fj^ is also zero in the quark target model of [29] . We 
also find that the contributions to f!^'"^^ from gluon fragmentation are power suppressed by Ph±/Q'^ or 
by as- This fact is in agreement with the factorization formulas in Eq.(9). In the below we will in turn 
show the factorization for F^^ and F^^ . 

3. We first study the factorization of f!^\ As mentioned in the above, the full calculation is very tedious. 
Instead of doing this we expand contributions with loop integrals in / first and then perform the loop 



F^\xB,Zh,Ph±) 
F!^\xB,Zh,Phi_) 



(fk±cfp± ( k± ■ 



cPk±cfp± ( kj_ 
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Figure 1: Diagrams for the contributions to the relevant structure functions. The broken hues are cuts. 
Black dots represent the insertion of the electromagnetic current. 



integral. A convenient way for the expansion is to analysis different regions of the loop momentum. The 
results should be the same as those obtained from the full calculation by taking the limit q\<^Q'^ . If one 
wants to have results of the whole g^-region to test the factorization, one has to do the full calculation. 

We denote the momentum of the gluon in the intermediate state as kg and the momentum of the 
virtual gluon as k. Each contribution from Fig.l can be written as: 



AL^^(2.'mi)(2.)H\P + , - p. - ■ ^' !■ ■ ■ 1 . (11) 



In each contribution there are five propagators, their denominators are denoted as Di[i = 1,2,- •• ,5). 
The numerator represented in the above as Tr [• • • ] is a trace of product of 7-matrices. We scale the 
momentum kg±_ = —Ph± as at order of A. and expand each contribution in A. Since the transverse 
momenta are at order of A, the leading order contributions comes from the integration region of k with 
k±_ at order of A or smaller, i.e., from the region with kj_ ~ (A, A) or k±_ ~ (A^, A^). Without performing 
the integration of k we can determine the order of A of each contribution. We use d to denote the leading 
power of a contribution by setting its numerator as 1 and n as the leading power of its numerator. 
Hence the leading order of the contribution is at A'^'*"'*. It seems that without a detailed calculation of 
numerators it is unclear how to determine n. But we know n > 1. It is rather easy to determine d of each 
contribution by power counting. We find that only those diagrams have the leading power with d = —4. 
They are the interferences of Fig. la to Fig. If with Fig.lg in two momentum regions. One is the Glauber 
region: 

F ~ (A2,a2,A,A), fc^ ~ (1,A^A,A), (12) 
where the virtual gluon is a Glauber gluon. Another is the ultra-soft region: 

/t^ ~ (a2,A2,A2,A2), /t^ ~ (1,a2,A,A), (13) 

By calculating each numerator of each contribution we find that the leading contribution comes only 
from the interference of Fig. la with Fig.lg with the virtual gluon as a Glauber gluon. It is at A~^. All 
other contributions to Fj. are at orders higher than A . It is constructive to discuss here what will be 
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factorized into Sivers function. Because of the outgoing gluon is collinear to the initial quark as indicated 
in Eq.(12), the Glauber region of the virtual gluon in Fig. la corresponds the collinear region of the gluon 
emitted by the initial quark. This gluon is collinear to the initial quark and has the final state interaction 
with the outgoing quark through an exchange of a Glauber gluon. We will find in our calculation of the 
Sivers function in Fig. 2a, the emission of the collinear gluon and the final state interaction is correctly 
contained in the Sivers function. It is straightforward to obtain the leading contribution: 



mal{Nlj-l) 5{l- Zh)xB{l- xb) 



In 



Svr^ \Phi.\ {PL + "^'(1 - ^b)') V^,\ + ^2(1 - xbY 



XB? 



(14) 



To verify the factorization formula in Eq.(9) we consider Sivers function of the same target as a quark 
with the momentum P and the transverse polarization s±. In order to generate nonzero /cj_ and an 
absorptive part, one needs exchange of two gluons at leading order. Nonzero absorptive parts in the 
amplitude can be found by cutting diagrams. The contributions at leading order of Us are given by the 
interference of diagrams given in Fig. 2. 



(a) 




(6) 



(d) 




(/) 




-h.c. 



Figure 2: Diagrams for the contributions to Sivers function. The broken lines are cuts. The double lines 
represent the gauge links. 

In Fig. 2 the diagrams in the first raw have an absorptive part indicated by the cut. It should be 
noted that the gauge link here is pointing to the future representing a particle in the final state with 
an infinity-large — component of its momentum. Hence the energy flow is from left to right. This is 
different than the case in Drell-Yan processes, where the gauge link represents an incoming particle with 
an infinity-large — component of its momentum. This difference results in that the cut diagrams for 
Sivers function in SIDIS are different than those for Sivers function appearing in Drell-Yan processes. 
This can be seen by comparing Fig. 2 here with Fig.l and Fig. 2 in [3j. Since we work at the leading order 
of as, we can take the limit u"*" ^ u~ directly here. After taking the limit, the Sivers function receives 
the nonzero contribution only from the interference of Fig. 2a with Fig.2g. This is also in correspondence 
to the case of The detailed calculation can be done similarly as done in [3]. We have: 

"5(^0 -1) mx{l-x) f m'^{l-xf \ 

q±{x,k^)- kUk^,+mm-x)^) (k^,+mHl ' ' ^ ' 



X 
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This result has been also obtained in [30]. Comparing the Sivers function for Drell-Yan processes in [3] 
we have the expected relation: 

q±{x, k±)\siDis = -q±{x, k±)\DY, (16) 

although it looks different from the corresponding cutting diagrams. The leading result of q{x,p±) can 
be found in 

q{x,P±) = 6{l-x)6\p^). (17) 

With the partonic F^^^ in Eq.(14) we verify the factorization formula for f!^^ at the leading but nontrivial 
order of Og- 

Before turning to F^^ , we briefly discuss the corresponding result for F^^ in the collinear factorization. 
The formal derivation can be found in [23\ I27j. where the perturbative coefficient in the factorization is 
obtained. With our partonic result for F!^^ in Eq.(14) and that for the twist-3 matrix element calculated 
in [3], we find that the collinear factorization for f!^^ is the same as for Drell-Yan processes given in [3], 
but with a perturbative coefficient different than that derived formally in [27j. In this factorization a 
twist-3 matrix element Tp{xi, X2, IJ-) defined in [20l [21] appears. With the result of Tf{xi, X2, fi) defined 
in [3] we have the collinear factorization: 

F^\xB,zh,Ph±) = 77^1^ / ——-D(y2)(^(l - 6)7; — ^-r^XBTF{xB,yi,Ph±), (18) 

27r^P^j_ J yi y2 (1 - 6)+ 

with ii = xs/yi, ^2 = Zh/y2- In the above D is the standard fragmentation function in the collinear 
factorization. The same situation is also found in Drell-Yan processes [3]. With the result presented in 
this work it is difficult to see why F^"^ derived from our partonic results takes a different factorization 
form than that derived with the diagram expansion at hadron level. In the collinear factorization for 
SSA [201 12H [22 I I23|. there are three partons entering hard scattering in cut diagrams for the differential 
cross section. We have constructed a suitable state of massless partons in our later work [34] to study 
the difference, where there are also three partons entering hard scattering. We have found that the 
so-called soft-pole contributions do not exist at the leading order of order and the gluon among the 
three partons is always transversely polarized. This can be the reason for the difference. Details can be 
found in [34] . 



4. Now we turn to Fj. . In this case it becomes complicated because more than one diagram can give 

nonzero contribution and also the loop integrals have a collinear divergence. Since our purpose is to verify 

(2) 

the factorization, we can first find the leading contributions to Fj. from Fig.l without performing the 
loop integrals, and then find the corresponding contributions to Collins function to check the factorization. 
The possible contributions to Collins function are given by the diagrams in Fig. 3. In Fig. 3 each of these 
diagrams is in one-to-one correspondence to each of those in Fig.l, if one identifies the gauge link in Fig. 3 
as the initial quark line in Fig.l. The result for 6q± at leading order simply reads: 

SqT{x,k^) =6{l-x)6'^{k±). (19) 

Let us look at the contribution from the interference of Fig. la with Fig.lh. We denote the momentum 
of the gluon crossing the cut in Fig. la as k. In this contribution the possible dominant contribution is 
with kg in the region collinear to Ph, i.e., kg is scaled as (A^, 0(1), A, A). With the power counting of 
the denominators and the integration measure of k we find the possible dominant contributions are from 
two regions of k. One is specified with k as collinear to kg, another is with k scaled as (A'^, A^, A^, A^). 
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Figure 3: Diagrams for the contributions to Collins function. The broken lines are cuts. The double lines 
represent the gauge links. 

Through calculating the corresponding trace in the numerator we find that the leading contribution to 

(2) 

W^'^ is only from the collinear region. We extract the structure function Fj^ as: 



Fr 



(2), 



I ah 



1 



1 



\P, 



h±\ 



W) 



P,k- {kg - kY + ie 



^6{k^)5{{Ph + kg-kf -m^) 



k- {2k; - k- + 4P, ) - p. 



+ • 



(20) 



where the • • • denote the terms suppressed by the power A or at higher orders of m . The leading order 

(2) I 

of Fj. \ah is at A with d = —4 and n = 1. In the above integral one can use the two (5-functions to 
perform the integration of A;^ and k~ . The integration over k± is bound from the above because of the 
momentum conservation. We can extend the upper bound to infinity because we only need the leading 
order contribution. We also note that the above integral has a collinear divergence when k is collinear 
to kg. The Collins function receives a contribution from the interference of Fig. 3a with Fig.3h. This 

(2) 

contribution corresponds to the above contribution Fj. \ah- Evaluating this contribution we have: 



{Zh,P} 



h±)\ah 



An 

1 



■(AT^-l)- 



d^k 



PL J (2-)2 



5{k'^)S{{Ph + kg - k)"^ - m"^) 



Prk- {kg - kf + ie 



k' {2k; -k~ + 4P, 



_ k± ■ kg± 

h ) 1.2 



{2kZ + k- 



(21) 



(2) 

Comparing the above equations we find that the two contributions satisfy the factorization of Fj. in 
Eq.(9). 

(2) 

Inspecting other contributions to Fj. ' in Fig.l, we find that the following contributions are at higher 
order of A: The interference of Fig.lg with Fig. la, Fig. Id, Fig.le and Fig. If and the interference of 
Fig.le with Fig.lh. Calculating the correspond diagrams in Fig. 3 in the limit ^ , we find that 



8 



the contributions to Collins function, corresponding to the above interferences, are zero. Calculating 
the remaining diagrams in Fig.l similar to the interference of Fig. la with Fig.lh, we find that the 
leading contributions can be factorized with the contributions to Collins function in Fig. 3 in an one-to- 

(2) 

one correspondence. Therefore, the factorization for Fj. in Eq.(9) is verified with our partonic state, 

(2) 

although we have not obtained detailed results for Fj, and the Collins function with our parton states. 

In e'*"e~ -annihilations into hadrons, certain SSA is related to the Collins function which is different 
than the Collins function in SIDIS. In e"*'e~ -annihilations the Collins function is defined by replacing the 
past-pointing gauge links in Eq.(8) with the future-pointing gauge links. It is interesting to ask if there 
is a relation between the two Collins function like the case with Sivers functions in Eq.(16). The relation 
for Sivers functions in Eq.(16) can be derived in general with Parity(P)- and Time(T)-reversal symmetry. 
However, it seems impossible to obtain a relation between two Collins functions with PT symmetries, 
because after PT-transformation the semi inclusive sum of out-states becomes the semi-inclusive sum of 
in-states. It is unclear how to relate a matrix element with such an in-state to that of the corresponding 
out-state[10j. In [lOj it is proposed that one can use the future-pointing gauge links instead of the past- 
pointing gauge links in fragmentation functions to perform the factorization. The proposal has been 
examined at one- loop level and arguments beyond one-loop level have also been given [10]. This implies 
that the direction of the gauge links is irrelevant in the two Collins functions and they are the same. 

Some model-calculations also have shown that the two Collins functions are the same[3H 1321 133j. In 
|33j it has been shown that the different directions of the gauge links in the case with one- and two gluon 
exchange make no difference in contributions to Collins functions. From our calculation we can also 
conclude at the considered order that the Collins functions are universal. The contributions in Fig. 3 to 
the Collins function in SIDIS are with the past-pointing gauge link. This gauge link represents a particle 
with infinite +-energy moving from right to left in Fig. 3. If we change the gauge link into the future- 
pointing gauge link, the future-pointing gauge link represents a particle with infinite -|— energy moving 
from left to right. With the future-pointing gauge link, the cut for Fig. 3b, Fig. 3c and Fig. 3d is the same 
and also there is no additional cut for these diagrams. Hence the contributions from Fig. 3b, Fig. 3c and 
Fig. 3d are the same for two Collins functions. The contribution of Fig.3e to Collins functions is always 
zero. With the future-pointing gauge link, the two remaining diagrams, i.e.. Fig. 3a and Fig.3f can have 
3 cuts. One of them is the same as given in Fig. 3. This cut gives the same contribution. The additional 
two cuts are cutting the gauge link and other propagator. It is interesting that the contributions from 
the additional two cuts cancel each other. Similar cancelation of contributions from additional cuts by 
changing the direction of gauge links is also observed in a model-calculation [31J . Therefore, our result 
shows that the two Collins functions are the same at the considered order. 

5. In this work we have verified the TMD factorization for SSA in SIDIS with partonic states at the 
leading order of Og- Although it is at the leading order, but it is nontrivial as shown in this work. Two 
structure functions related to SSA has been examined. They take a factorized form in terms of TMD 
parton distributions and TMD parton fragmentation functions. The factorized form agrees with that 
derived formally. However, as in the case of Drell-Yan processes, we find that our partonic results satisfy 
a collinear factorization but with a perturbative coefficient different than that derived formally. In this 
work and [3] we have made attempts to check two factorization approaches with partonic states with a 
nonzero quark mass. It is worth to point out that the same results of SSA related to Sivers function can 
also be derived with massless partons as shown in our later work [3^. 
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